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Enhanced Dust Emission in the HL Tau Disc: A Low-Mass 
Companion in Formation? 
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ABSTRACT 

We have imaged the disc of the young star HL Tau using the VLA at 1.3 cm, 
with 0.08" resolution (as small as the orbit of Jupiter). The disc is around half the 
stellar mass, assuming a canonical gas-mass conversion from the measured mass in 
large dust grains. A simulation shows that such discs are gravitationally unstable, and 
can fragment at radii of a few tens of AU to form planets. The VLA i mage shows a 
compa ct feature in the disc at 65 AU radius (confirming the 'nebulosity' of lWelch et alj 
(|2004l) ). which is interpreted as a localised surface density enhancement representing a 
candidate proto-planet in its earliest accretion phase. If correct, this is the first image 
of a low-mass companion object seen together with the parent disc material out of 
which it is forming. The object has an inferred gas plus dust mass of s=s 14 Mj up ; ter , 
similar to the mass of a proto-planet formed in the simulation. The disc instability 
may have been enhanced by a stellar fly by: the proper motion of the nearby star XZ 
Tau shows it could have recently passed the HL Tau disc as close as ~ 600 AU. 

Key words: planetary systems: formation - planetary systems: protoplanetary discs 
- stars: pre-main-sequence - circumstellar matter - radio continuum: stars. 



1 INTRODUCTION 

The mechanisms by whi ch giant planets form are uncertain. 
Core - accretion models (|Pollack et al. I ll996l : iHubickvi et al.l 
[2005j, e.g.) have successfully linked high abundances of 
rocky elements in the star to higher pl anet-probability 
(|Fischer fc Valentil 120051 ; ISantos et al.ll2004 e.g. ), and may 
expla in planets with substantial rocky cores (iSato et al.l 
120061 ) - but have theoretical difficulties with slow planetary 
cooling that limits mass accretion rates and with rapid in- 
wards migration leading to loss of cores into the star. Also, 
the time of around 6 Myr to complete Jupiter may conflict 
with the infrared detection rate of dis cs that declines close 
to zero by 6 Myr (|Haisch et al.ll200lD . and with the latest 
ages of ~ 15 Myr when gas is detected (|Dent et al J 120051 ) 
since planet completion takes longer in low-mass discs. The 
alternative model of gravitational instability can create a 
proto-planet very rapidly, on dynamical (orbital) timescales 
i|Bosslll997l ; iRice et al.ll2005h. provided that the disc cooling 
time is similarly short i|Gammidl200ll ; lRafikovl l2005'). Only 
relatively rare discs of > 0.1 x M s tar will be unstable, but 
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recent radio studies (|Rodmann et al . 2006], e.g) that account 
for mass in large dust grains may boost more of the total 
disc masses into this category. Since me tre-sized particle s 
are gathered up into the gas fragments l|Rice et all 120061 ). 
this model may also account for solid cores to giant planets. 

Imaging a planet within its birth-disc would illumi- 
nate the processes involved. Companions of several Jupiter 
masses upwards have been imaged, but large separations 
from the primary and in some cases small mass-ratios of 
the two components sugg est that these objects may have 
formed like binary stars (jLuhmann et all |2006| ). No discs 
within these systems have been imaged, so the formation 
processes remain obscure - the most direct observational 
evidence for a fo rming object is a cleared cavity within the 
disc of AB Aur (|Oppenheimer et al.ll2008l ). 

In this letter, we present the first candidate for a low 
mass companion imaged in the accretion stage and within 
the parent disc. Such an object would be expected to ap- 
pear as a cool condensation of dust and gas, possibly with- 
out a distinct dense core as sediment ation timescales fo r 
large dust grains can be a few 10 4 years (|Helled et alll2008h . 
Here, we use radio- wavelength data to trace the thermal 
emission from large dust particles in the disc around HL 
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Figure 1. VLA 1.3 cm images towards HL Tau. Main image: natural weighting with a beam (small inset) of 0.11" for maximum 
sensitivity; contours increase by \f2 and are at (4.0, 5.7, 8.0, 11.4, 16.0) xl7|xjy/beam (lcr). The ar row indicates the je t axes and the 
ellipse shows the approximate extent of the inclined disc in a previous 0.6" resolution image at 2.7 mm l lLoonev et al . 2000). The compact 
object lies to the upper-right. Upper inset: same but with central peak of 263 ujy subtracted, to highlight the jet bases and two features 
at fa 20 AU at the ends of the disc major axis. Contours from the unsubtracted image are overlaid. Lower inset: uniform weighted image 
at higher resolution of 0.08", with contours at (3.0, 4.2, 6.0, 8.5, 12.0) x21ujy/beam (lcr), highlighting the compact object. 
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Tau. This pre-main-seque nce Class I (remnant en velope) ob- 
ject has been modelled bv lRobitaille et al.l \200H ) at around 
0.33 M Q and 5 L Q , seen at < 10 5 years old. The HL Tau 
disc was selected as one of the brightest known at mil- 
limetre wavelength s, with estimates for gas plus dust mass 
of up to 0.1 M Q l|Beckwith et al.lll990l ). and thus within 
the disc-to-star mass regime where instability could occur. 
Millimetre interferometry (resolving out the envelope) has 
shown emission fro m the dust-disc extending out to a t least 
~ 100 ALT rad i us (Iwilner et al.1 Il996l; iMundv et al.|[l99r3 ; 
lLav et al.lll997l ; lLoonev et alj|2000l : iRodmann et al-lbOOrf T 



2 OBSERVATIONS 

Observations of HL Tau in a 10" field centred at RA 
04:31:38.4034, Dec. 18:13:57.748 (J2000) were made with 
the Very Large Array at 1.3 cm wavelength. In the largest 
A-configuration plus the 50 km-distant Pie Town antenna, 
the VLA was sensitive to scales down to 0.08", a factor 
of three higher tha n the best previous resolution of 0.25" 
l| Welch et al.ll2004l ). At the Taurus distance of « 140 pc, 
this gives a resolution of just over 10 AU, equivalent to 
the orbit of Jupiter. The 22.5 GHz data were obtained 
over two runs in 2006 Mar and Apr, with a total usable 
time of ~ 12 hours. We used the phase reference source 
04311+20376, 2.4° from HL Tau and the primary flux scale 



was provided by 3C286. We also corrected the antenna 
pointing and refined the amplitude calibration with the aid 
of bright compact sources including 0552+398. We followed 
standard VLA observational and data reduction procedures 
as described at http://www.vla.nrao.edu/astro/ including 
procedures for High Frequency Data Reduction. Natural 
weighting gave a beam size of 114 by 104 mas and la sensi- 
tivity of 17 ujy per beam. The data were also reconstructed 
with uniform weighting, giving a smaller beam of 82 by 76 
mas with la sensitivity of 21 ixjy/beam. Systematic posi- 
tional uncertainties are ~ 30 mas, less than the resolution. 

We also observed HL Tau at 5 cm wavelength with the 
MERLIN array (using up to 6 antennas) in 2006 Jan-Feb 
for a usable total time of 20 hr including calibration. We 
used the phase reference source B042 5+174 at 1° separatio n 
and followed procedures described in lDiamond et aL l (|2003h . 
These were some of the first observations made with six 5- 
GHz receivers (not all cryogenic), reaching a la sensitivity 
of 100 ^iJy beam -1 using a 100- mas restoring beam. The im- 
ages were sensitive to scales of 0.04-0.8 arcsec depending on 
weighting. We re-observed in 2007 Apr at 6-cm wavelength; 
the combined data reached a sensitivity of 55 ujy beam - . 
This gives a formal 3a upper limit for 5-6-cm emission on 
these scales of 165 ujy beam -1 . No emission brighter than 
100 fiJy beam -1 was detected within 0.2 arcsec of the com- 
pact object discussed in Section 3. 
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3 RESULTS 

The image (Figure 1) shows a bright inclined disc out to 
around 30 AU, i.e. similar to Neptune's orbit. The elliptical 
morphology and similar orientation to lower-resolution im- 
ages confirm that this is emission from the dust disc, while 
the two orthogonal extensions (upper inset) are presumably 
the base s of the bipolar jet, previously seen at 1 -2" from 
the star (|Rodmann et alj 2006; Lo onev et alj|200(il ). Fainter 
dust emission is detected out to 100 AU, declining with ra- 
dius as r _1 - 5 ± - 5 (from an error-weighted fit using elliptical 
annuli from 0.1 — 0.7", with a correlation coefficient of 0.9). 
Assuming emission weighted by r _1//2 for grains in thermal 
equilibrium with the star, the disc surface density decl ines at 
r~ 1± 0,5 , roughly like the young outer Solar System l|Davisl 
120051 ). Imaged features include the extended disc (~ 1350 
u.Jy); a central peak of ~ 300 U-Jy located at (J2000) 04 
31 38.4184, +18 13 57.387 (2 mas fit errors); the exten- 
sions to the NE and SW interpreted as jets, each of flux 
ft* 100 ± 25 nJy; and a clump of 78 ± 17 u.Jy offset from 
the central peak by 380 mas at position angle (PA) of 30° 
clockwise. The jet, disc and clump features are highlighted 
in a higher-resolution image (lower inset). 

This compact clump is at a projected stellar separa- 
tion of 55 AU, or orbiting at around 65 AU if corrected for 
proje ction assuming a disc inclination of 60° l|Wilner et al.l 
1 19961 ). The clump is here resolved for the first time. A 'neb- 
ulosity' wa s previously repor ted in a BIMA 1.4 mm image 
of the disc (|Welch et al.ll2004h . separated by 70 AU from the 
star at PA -40°. Within the BIMA resolution of 0.25" (35 
AU), this is coincident with our VLA 1.3 cm peak at 55 AU, 
—30°. The two independent detections give high confidence 
that this feature is real, while it is seen here in detail for 
the first time. The earlier image of an unresolved flux en- 
hancement could have been attributed to an ordinary disc 
asymmetry such as a large-scale perturbation, but in our 
new data a clump is clearly seen. It is compact (full-width 
half-peak sizes of 20 ± 12 AU by ^ 12 AU), three times 
brighter than the local flux level of the disc, and clearly sep- 
arated (by five resolution elements) from the stellar position. 
Given the similarity to proto-planets formed in simulations 
(Figure 2), we propose this object as a candidate for the 
earliest stage of growth of a low-mass companion. 



3.1 Robustness 

We tested whether this +4.5<r feature could be an arte- 
fact. The noise/pixel has a very Gaussian distribution across 
the image, so the probability of a random fluctuation of 
^ 78(i.Jy/beam occurring within the projected disc area out 
to 100 AU should be only ~ 0.03 %. As a test, we added 33 
fake 'planets' to the visibility data, with the flux density of 
the observed clump, scattered at 1 — 3" offsets, and followed 
the standard imaging process. The recovered objects follow 
a roughly Gaussian distribution, with mean and dispersion 
in flux density of 75 ± 25 (i.Jy. This distribution was com- 
pared to the fluxes of 33 random positions, away from the 
disc, and no noise pixels at these positions exceed 3a. The 
faintest fake 'planets', in the lower la tail of their flux dis- 
tribution, were close to the off-source 3a noise level, showing 
that a true planet of such low intensity might be lost in the 



Figure 2. Example image from an SPH simulation (see text) 
showing the surface density structure of a 0.3 Mq disc around a 
0.5 Mq star. A single dense clump has formed in the disc (upper 
right), at a radius of 75 AU and with a mass of 8 Mj up ; ter . 

noise - but no false planet is seen emerging as a 3a artefact, 
let alone at the 4.5cr level of the actual feature in the disc. 

We then subtracted out the bright central emission 
(Figure 1, upper inset) and found that the compact object 
is unaffected; conversely we added a model of this emission 
at random positions across the field - no spurious features 
above 3a were seen at the distance of our candidate. Hence, 
the clump is not an artefact arising from imperfect deconvo- 
lution of the interferometric image. Finally we split the data 
into two separate frequency bands and alternatively into left 
and right circular polarizations. The clump was always re- 
covered at a similar intensity, at between 2.5<r and 4.8a sig- 
nificance depending on the quality in the partial dataset, 
supporting a real detection. 

We also investigated whether an unrelated source could 
be seen through the HL Tau disc - such background ra- 
dio objects would typically be active galactic nuclei. This 
is improbable as even the most distant extr agalactic faint 
radio sources are at least ~ 0.4" in diameter (|Muxlow et al.l 
2005), five times larger than our beam. Conversely, a source 
would need to be well over 1 mjy in lower-resolution sur- 
veys to be detectable, and as AGN emission rises at long 
wavelengths (opposite to dust), this would be a rare bright 
object. Extrapolating a 1 mjy 1.3-cm source with a spec- 
tral index of a < — 0.2 (for 92 % of r adio sources with 
multi-frequency data (|Volmer et al.|[2005h and flux <x v~ a ) 
yields a signal > 2 mjy at 20 cm. The VLA FIRST survey 
(http:/ /sundog. stsci.edu/first/ ) shows such an object would 
turn up in our 10" x 10" field in about 1 in 3000 cases 
(actually much lower since in this flux range most sources 
are 3> 0.4" in size). Moreover, such a source would have 
probably appeared at 5-cm: for a < —0.2 the counterpart 
to the 1.3 cm feature would be of > 100 ^Jy. The MER- 
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Figure 3. Spectral energy distribution of clump (squares) and 
integrated flux around HL Tau, includi ng VLA (cross), MER- 
LIN upper limit (arrow) and results from[Rodmann ct al. (2006) 
(filled ovals). Unfilled ovals show the minimum flux-contribution 
from dust, after subtracting free-free flux extrapolated from the 
MERLIN upper limit assuming F v oc u 1 , the steepest slope typi- 
cally seen for ionised stellar winds jAnglada ct al. 1998). A flatter 
{y~ 01 ) subtraction gives an improved v 2 5 dust power-law fit. 
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Figure 4. Plot of cooling time as a function of radius for a disc 
mass of 0.13 Mq around a star of mass 0.33 Mq. The disc sur- 
face density declines as r — 15 to an outer radius of 100 AU, and 
the instability parameter is Q ~ 1. The quantity t coo \Q is a di- 
mensionless quantity relating the cooling time to the local orbital 
period; for t coo if2 = 2n the cooling time equals the orbital period. 
The dotted lines show critical values below which frag mentation 
occurs for different equations of state iRice et alj|2005h . 



LIN image showed no features above 100 uJy in this region 
(3cr limit of 165 ujy), so a background synchrotron source 
is ruled out with high confidence. A millimetre counterpart 
(see below) indicates a dust-like spectrum for the clump - a 
distant starburst galaxy could potentially have such a spec- 
tral energy distribution, but these are rare. Extrapolating 
from the 1.3 cm flux would lead to a source of > 20 mJy 
at 0.85 mm wavelength, wit h a corresponding probability of 
< 10 -6 l|Coppin et al.ll2006D within 100 AU of HL Tau. 



3.2 Measurements 

The spectral energy distribution confirms that circumstellar 
dust is detected. The integrated centimetre flux was previ- 
ously thought to be from the ionised stellar wind, but Fig- 
ure 3 shows dust emission to A 3.6 cm. The dust spec- 
trum of F„ cx V 2 - 5 - 2 6 is characteristic of emission from a 
population of particles extending i n size up to a t least three 
times the observing wavelengths (|Drainell2006l ). and hence 
here to bodies of > 10 cm. The spectrum of the conden- 
sation is tx v 2 ' 5 from the fluxes of 23 ± 5 mJy at 1.4 mm 
and 78 ± 17 ujy at 1.3 cm (n eglecting lower-resolution data 
noted by IWelch et all {2004) as surrounding disc flux may 
be included), again implying that very la rge particles are 
present. This would agree with simulations (|Rice et al.ll2006l ) 
in which 'boulder'-like bodies of around metre-size are most 
readily captured in unstable regions. 

At 1.3 cm, the measured disc flux excluding the jets and 
region inside 5 AU is ~ 1350 ujy. Assuming simplistically 
that the dust particles are in ther mal equilibrium with th e 
star, we adopt the best-fit 5 Lq from [Robitaille et al.l l|2007l ). 
giving 90 K at a characteristic disc radius of ~ 20 AU (Fig- 
ure 1, upper inset) and 50 K at the clump orbit. For an 
opacity of 7xl0 -4 cm 2 /g a t 1.3 cm (for populations extend- 
ing up to 1-10 cm particles (|Drainel2 006). and assuming that 
all the original gas and dust are present so that a canoni- 
cal gas-to-dust mass ratio of 100 applies), the disc contains 
around 0.13 Mq in total. This exceeds previous estimates 



of up to 0.1 Mq, due to the additional large dust and con- 
sequent scaling up of the total mass. The clump comprises 
~ 14 Mjupiter, with the uncertainty dominated by the scal- 
ing with temperature; the errors in flux and distance con- 
tribute at up to ~ 30 % levels while adopting a different 
opacity (e.g. v 5 extrapolation from a standard 0.01 cm 2 /g 
at 1 mm) reduces the mass by factor of two. The central 
peak may also contain a mass reservoir - the 1.3 to 6 cm 
spectrum is just consistent with a wind origin, but allows 
up to ~ 13 Mj U pit cr of gas and dust to be present if the 
1.3 cm flux is dust-dominated. This is si milar to th e pri- 
mordial 12 Mj U pitcr out to Jupiter's orbit (|Davis!l2005h . and 
so planets might form here by core-accretion; if large grains 
are present this may also solve the planetar y radiative cool- 
ing p roblem by reducing the dust opacity dHubickvi et al.l 

3.3 Simulations 

We ran an example 250,000-particle Smoothed Particle Hy- 
drodynamics (SPH) simulation to investigate if the HL 
Tau disc could be gravitationally unstable (Figure 2). The 
disc is assumed to have an initial surface density pro- 
file of E tx r -1,5 and initial outer radius of 100 AU. 
The simulation evolves under a radiative transfer formalism 
dStamatellos et a 1 2007) in which each particle is allowed 
to cool towards its equilibrium temperature, determined us- 
ing the local optical depth (dependent on material opacity, 
here taken to be half- Solar). The disc is allowed to heat up 
through pdV work and viscous dissipation. Relative to the 
local dynamical timescales, the cooling times are slow in the 
inner disc where the optical depth is high and faster in the 
outer disc. Fragmentation is expected if the cooling time is 
less than a fe w orbital period s, and depends on the equa- 
tion of state |Rice et al.| [2005). For the estimated star and 
disc masses of HL Tau, Figure 4 shows that the cooling time 
should be fast enough for fragmentation beyond ~ 40 AU. In 
the simulation, we adopted somewhat higher star and disc 
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masses (but within the uncertainties) and allowed the disc 
to cool as low as 10 K, which favours fragmentation (Fig- 
ure 2). A single dense clump has formed in the disc, and is 
located at a radius of 75 AU and has a mass of ~ 8 Mj up iter, 
although it may continue to accrete from the disc. These 
properties are similar to those of the actual observed object 
around HL Tau. The simulation results will vary with the 
chosen opacity - a smaller value produces additional clumps, 
while a larger one could inhibit fragmentation altogether - 
and with the surface density profile - a flatter profile yields 
more outer-disc mass and higher tendency to fragment. 



4 DISCUSSION 

Including the large grains now detec t ed, th e HL Tau disc 
mass is ~ 0.13 Mq. iRobitaille et~ai1 (|2007l ) find good fits 
to the stellar mass for 0.2-1 M©; for their best-fit value of 
0.33 Mq the disc is around 0.4M star . This proportionally 
massive disc should be gravitationally unstable, and a simu- 
lation at the higher end of the Mdi sc ,star ranges confirms that 
planetary objects could form at a few tens of AU. The VLA 
data show such a flux peak in the parent disc material, inter- 
preted here as a surface density enhancement. (The clump is 
three times brighter than the local disc flux, while warming 
of the gas by gravitational collapse should only contribute 
marginally to higher emission; simulation results suggest the 
beam-averaged dust temperature is raised by « 50 56.) This 
clump at 65 AU from HL Tau lies in the appropriate unsta- 
ble region, and is compact as expected for a low-mass object 
accreting from the disc. 

The simulated disc is unstable for the adopted parame- 
ters, but external forces could have increased the real disc's 
tendency to fragment. Notably, another cluster member, 
XZ Tau, appears close-by, which is unusual within the dif- 
fuse Taurus association, and the relative motions suggest 
a possible recent encounter of the two stars. Their line-of- 
sight distances are unkn own, but the similar radial velocities 
(Fol ha fe Emersonl2 000') suggest they are not located in very 
different parts of the association, and in 2-D the stars are 
presently divergi ng. XZ Tau lies 23" eas t of HL Tau and the 
proper motions l|Ducourant et ai] |2005h are (+11,-19) and 
(-3,-21) milliarcsec/year respectively (errors of 2-5 mas/yr). 
Around 1600 years earlier, the stars could thus have passed 
within ~ 600 AU (in 2-D projection). Such an event would 
have been dynamically recent, given that the compact object 
has an orbital period of 900 years for M s tar of 0.33 Mq. 

The final mass of this still-forming companion may in- 
crease, by absorbing more of the disc, but our estimate of 
14 Mj U pitcr for the condensation is well down into the sub- 
stellar regime. If all this material is accreted, the final object 
would be around the brown dwarf / planet boundary by the 
definition of short-lived deuterium-burning capability, which 
occurs at > 12-13 Mj up iter. A more recently developed defi- 
nition of a planet is a low-mass object that formed in the disc 
of a star. This 'origins' definition sides t eps th e deuterium- 
burning issue, which as IChabrier et ail (|2007f ) point out is 
irrelevant for the evolution of brown dwarfs. In the case of 
HL Tau 'b', imaging the object within the parent disc marks 
it as a candiate proto-planet by this origins definition. 
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